Stable Elliptical Vortices in a Cylindrical Geometry 
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We show that, in a two-dimensional (2d) ideal fluid (also applies to a column of quasi-2d non-neutral 
plasma in an axial magnetic field), large elliptical vortices in a finite disk are stable. The stability is 
established by comparison between energy of elliptical and symmetrical states to satisfy a sufficient 
condition, without dynamical eigen-analysis. Analytical small ellipticity expansion of energy and 
exact numerical values for finite ellipticity are both obtained. The expansion indicates stable linear 
I = 2 diocotron modes for large vortices (or plasma columns). Numerical simulations of the 2d Eulcr 
equation are also performed. They not only confirm the sufficient condition, but also show that the 
stability persists to smaller vortex sizes. The reason why decaying I — 2 modes were obtained by 
Briggs, Daugherty, and Levy [Phys. Fluids 13, 421 (1970)] using eigen-analysis is also discussed. 



The two-dimensional (2d) incompressible Euler equa- 
tion 



da; 

— + (u • V)w = 0, 



(1) 



not only describes an incompressible 2d ideal fluid, but 
also governs the behavior of a long non-neutral plasma 
column confined by a uniform axial magnetic field |Q. 
Here u(x, y) is the 2d velocity field and u>(r) is the vortic- 
ity field, uj = (V x u) -z. The incompressibility condition, 
V • u = 0, can be automatically satisfied by defining the 
stream function <f> as u = (d<p/dy, —d<fi/dx). The stream 
function and vorticity are related by the Poisson equation 
V 2 = — uj. In a pure electron plasma, ui corresponds to 
the electron density and <f> to the electrical potential. 

Stability problems of coherent vortex states in this 
system are long being interesting and important ques- 
tions. In a free space, there exist exact nonlinear el- 
liptical (Kirchoff) patch solutions g]. In a cylindrical 
geometry Briggs, Daugherty, and Levy [|| showed that, 
using dynamical eigen-analysis, resonance between fluid 
elements and wave modes will lead to damping of Z > 2 
diocotron modes. Here I denotes the mode number as 
the perturbation to a symmetric stream function is writ- 
ten as 4>i(r) exp[z(f2t — 18)]. By solving the initial value 
problem of linearized equations and properly treating an- 
alytical continuation in complex Q plane, they obtained 
formulation for complex eigenvalue f2. In particular, for 
a vorticity distribution very close to a step function but 
negative radial derivative at all places, f2 with a positive 
imaginary part is calculated for I > 2, leading to decaying 
normal modes. 

Experimental observations of decaying I = 2 modes 
have been performed by Pillai and Gould || in a pure 
electron plasma. Exponential decay rates were obtained, 
as well as the observation of fluid trapping in the dio- 
cotron mode at large amplitudes. In another experiment 
with a pure electron plasma, || beat-wave resonance 
damping (transitions from high I modes to low I modes) 
was observed to be the dominant vortex symmetrization 
mechanism. 



A stability argument based on global constraints has 
also been applied to the 2d vortex system The logic 
of this analysis is to show that a functional W[uj] which 
is conserved by the 2d Euler equation is a maximum at 
a particular uj(y) against all other states that are ac- 
cessible under incompressible flows. At this maximum, 
no further changes in uj(r) are possible and the state is 
then stable. For example, Davidson and Lund ]TJ showed 
that a state in a cylindrical geometry following a relation 
oj(r) = u>((f>(r)) and du>{(j>) / d<j> > is nonlinearly stable 
||. In another example, O'Neil and Smith |J demon- 
strated that an off-center coherent vortex (linearly an 
I = 1 perturbation) in a disk is also stable. However, 
no results on the stability of an I — 2 mode using this 
method have been given in the literatures. 

Thermal equilibrium has been studied in 2d ideal fluids 
|l0| . Since the coarse-grained entropy will not decrease 
due to the dynamical vorticity mixing, it is proposed that 
the system will reach a maximum coarse-grained entropy 
state at long time. Mean field equations governing these 
states have been derived JToj] , and solutions in some situ- 
ations were obtained pd[ | . Once a mean-field equilibrium 
state is obtained, its stability can be assured by showing 
a positive second derivative of entropy against all possible 
perturbations. This of course is using a similar principle 
as the method mentioned above. 

In this paper we establish the stability of a large ellip- 
tical vortex (comparing to the system size) against relax- 
ation to a symmetrical state using neither of the above 
two methods with eigen-analysis and global maximum. 
We will first deduce a stable sufficient condition and then 
show that it is satisfied by elliptical vortices larger than 
a critical radius. The method is to compare energy of 
proper states, not by evaluating second derivatives, actu- 
ally not even finding any equilibrium states. We further 
perform numerical simulations of the 2d Euler equation 
to test our predictions. Simulations not only confirm the 
sufficient condition, but also show that elliptical vortices 
are stable to lower radii. 

Basic argument of the sufficient condition goes as: 
Consider initially a uniform- vorticity elliptical vortex sit- 
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ting at the center of a unit disk, with unit vorticity level 
without losing generosity. Now consider its possible dy- 
namics toward an axis-symmetrical vortex. This will be 
a state with a linear I = 2 diocotron mode if infinitesimal 
cllipticity. 

The Euler equation conserves the total vorticity Q, an- 
gular momentum M, and energy E of the initial ellipse, 
which are given by 

Q = J uj(r)dr M = J r 2 uj(r)dr 

E e = \ f <t>(r)w(r)dr 

Furthermore, dynamical vorticity mixing ensures that 
the vorticity level of the resulting symmetrical vortex will 
never exceed one (the original uniform value). Under this 
restriction and given Q and M from the initial ellipse, 
there must be a maximum energy state with its energy 
denoted as E s among all possible symmetrical distribu- 
tions. With conservation of energy, this condition then 
immediately follows: 

E e < E s is necessary for the ellipse to ever 
evolve to a symmetrical vortex; 

E e > E s is the sufficient condition for the 
ellipse not evolving to a symmetrical state. 

Applied to infinitesimal ellipticity, the I = 2 diocotron 
mode will not decay when E e > E s . 

It should be noted here that this condition only try 
to exclude symmetrical states from possible evolutions, 
a limitation purely physically motivated. For example, 
it seems unlikely that an ellipse at the disk center will 
break the symmetry and relax to an off-center vortex, 
although we believe that the energy of off-center vortices 
could be larger than E e . The conjecture (not decaying to 
off-center states) is confirmed by numerical simulations 
which will be discussed later. 

To test the above condition, our first task is to calcu- 
late the energy E e of a uniform elliptical vortex, which 
we define as a vorticity distribution Lu e (r) in the polar 
coordinate (r, 9) 



Lu e (r,9;r 0l e) = 1 - s(r - r (l + e cos 26*)), 



(2) 
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We separate E e into three terms in the last equation. 
Here <f>o(r; ro) is the stream function of a uniform circular 
vortex with radius r , ^(r^-) = -wo, w (r;r ) = 
1 — s(r — r ) and E its corresponding energy, 



E n = 



(f>o(r;r Q )uJo(r;ro)2irrdr = Trr 4 , (-±mr + ^) 



We know of no way to integrate Eq. (^) analytically. 
Nevertheless, we can study the linear stability of an / = 2 
diocotron mode from small e behavior of E e (rQ, e). Since 
the vortex is defined by ro(l + e cos 28), the lowest order 
dependence on e must be e 2 . Correct to the order of e 2 , 
the second term in Eq. (0) is quickly found to be, 
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I [ro</>o(r ;r ) + M r o] r )} r%e 2 cos 2 26d6 



= -i7r^(l + lnr )e 2 . 

Here prime denotes the derivative respected to r. Evalu- 
ation of the third term in Eq. (^ is more difficult. Again 
correct to the order of e 2 , the integration becomes 
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2tt 
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d9' cos 29 cos 29'G(r ,9;r ,9'). 



Using the Green function and changing to new variables 
u = 9 + 9',v = 9 — 9', after some algebra, we reach 
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ln(a — cos v) cos 2vdv 



(4) 



with s(x) the usual step function. The parameter ro de- 
fines a base vortex size and e its ellipticity. To compute 
the energy of this vortex in a unit disk, first we consider 
the Green function in a disk for the Poisson equation, 
V 2 G(r;r') — —S(r — r'), with zero boundary condition 
at r = 1. Using an opposite-charged image charge sitting 
at r" = (l/r',9 1 ), the Green function can be written as 
G(r;r') = (In |r - r'| -ln|r-r"| -lnr'). The en- 
ergy of the uniform elliptical vortex is then 



E e (r ,e) 



(r)tu(r)dr 



with a = Tj(r 2 + 1/r 2 ) > 1. The integration / = 

J 27r ln(l — cos u) cos 2udu = — n has also been used in 
reaching Eq. (|J). 

The integration in Eq. (Q) is computed by first inte- 
grating its derivative respective to a, and then using / 
to determine the constant arising from integration of a. 
Eventually Eq. (||) is found to be |7rr 4 (l — r 4 )e 2 , and 
the energy of the elliptical vortex becomes 



E e (r {) ,e) = E + {nr 4 Q ( 



-±r 4 

2'0 



lnr ) e 2 + 0(e 4 ). 



(5) 
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The energy E e (ro,e) is now to be compared with the 
energy E s of the maximum-energy symmetrical state 
with the same values of total vorticity Q and angular 
momentum M. Its vorticity must also be equal or less 
than one. To see what this state is, first it is favorable to 
have all the vorticity stay together, i.e., a uniform unit- 
valued circular vortex with radius r s — (Q/n) 1 / 2 , to gain 
as much as energy. However this circular vortex has a 
fixed angular momentum \t^t\ , and the uniform ellipse 
always has a larger value. To satisfy the requirement of 
both Q and M, as well as achieving a maximum energy 
the state will have a vorticity distribution u) s (r) as, 

, , fl for 0<r<a & /?<r<l 
^M = ( for a<r<f3. (6) 

Here a and (3 depend on Q and M, which are determined 
by ro and e. In this distribution, a certain amount of 
vorticity is put as far away from center as possible, i.e., 
at the disk boundary, to account for the excess angular 
momentum and maximum amount of vorticity is left to 
concentrate at the center to acquire a maximum energy 
Jl2| |. Here we see how the system size comes into play 

in a delicate manner. At small e, a = (l + | e 2 ^) 

4 

and (3=1 — tjt-^tC 2 ', and the energy E s is expanded as 
(It involves only straightforward algebra to solve 4> s and 
then integrate E s .) 

E s =E,-\^ l —^l-\nr^ + 0{e i ). 
4 1 — ?'g 

Now we obtain the energy difference between E e and E s 
as 

Ee - E s = \rt (-\rt - I - ^ lnr ) e 2 + 0(e% 

Evaluation of e 2 term reveals that there is a critical value 
of ro, r c w 0.586, such that E e < E s for ro < r c and 
E e > E s for r > r c . 

So applying the energy condition, this indicates: the 
I = 2 mode perturbation of a circular vortex in a finite 
disk will not decay if the vortex is large enough (larger 
than 0.586 times the disk radius). This result seems con- 
tradict that of Briggs, Daugherty, and Levy || where 
decaying modes were calculated from eigen-analysis for 
alH > 2 modes of a circular vortex with a smooth profile 
very close to wo(r;ro) (a step at ro) but negative u)'(r) 
at all r. The resolution lies at that in the calculation of 
Briggs, Daugherty, and Levy, the symmetrical vortex is 
assumed as a monotonic decreasing function of r. This 
seems a reasonable and harmless condition. However, 
as Eq. (g) shows, this condition is very restrictive and 
always violated by the uniform ellipse and hence their 
results no longer apply. 

To further determine the stability of an ellipse with 
finite ellipticity, we need go beyond the expansion and 
calculate the energy for arbitrary e. Here we resort 



to numerical calculations of integration E = ^ J (pujdr. 
Since the Green function using image charges has loga- 
rithmic functions and is not easy to handle numerically, 
we rewrite the Green function as a summation of Fourier 
components in the azimuthal direction, 

oo 

G(r; r') = £ g m (r, r') cos(m(0 - 9')), 

m— 

with g m power ±m of both r and r' . The energy now be- 
comes a summation on m of four-dimensional (r, (9, r' , 8') 
integrals. The integration on r and r' can be carried out 
analytically and the energy simplifies to a summation of 
double integrals on 6 and 6' . The integrals are then cal- 
culated numerically, and results are checked to conform 
to Eq. (j|) at small e. 

The exact value of E e (ro, e) now enable us to establish 
the stability of finite ellipticity. In Figure 1 of the ro~e 
plane, we plot a solid line indicating the position where 
E e = E s . (Although complicated, again E s with arbi- 
trary e can be written down analytically from straight- 
forward algebra). To the right of the line, E e > E s and 
an elliptical vortex will never relax to a symmetrical vor- 
tex. The line is almost vertical and only curves a little to 
the left as e is increased. It crosses e = at ro « 0.586, 
the value we have obtained from the small e expansion. 
To the left, the present analysis only says that the decay 
to a symmetric state is allowed, but its occurrence is not 
implied. 
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FIG. 1. The vortex size and ellipticity space. The solid 
line marks the position where E e = E s . Squares represent 
relaxations to elliptical states in simulations, and circles to 
symmetrical states. 

It should be emphasized here that we have proved that 
the ellipse defined by Eq. (||) will not decay to a symmet- 
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rical state if r$ > r c . It is very likely that dynamically 
it will undergo adjustment and reach an elliptical-like 
steady state. We cannot say about its exact distribu- 
tion. Current understanding is that it probably should 
be a state described by u(r) = u>(4>(r) + fir 2 ), with 20 
giving the rigid body rotation frequency around the disk 
center. With a particular assumption on this functional 
dependence, an exact distribution can then be computed. 
One example is the mean field equilibrium JlO[ . However, 
whether and when the system will reach the prediction 
from this maximum-entropy principle (thermal equilib- 
rium) is still not very clear |l~3| ]. 



Although Eq 
pect that a smoot 



defines a uniform vortex, we do ex- 
;hly distributed vortex should have sim- 
ilar stability property if not deviated too much from Eq. 
(||). The exact values of stable radii will of course be 
different. This is confirmed by the numerical simulations 
discussed next. 

We also perform numerical simulations to test our pre- 
dictions. Simulations of the Euler equation in the polar 
coordinate have the difficulty of singularity at the origin 
due to vanishing grid spacing. To avoid this singularity, 
we use the functions, 



a; = /Vl-C72 

mapping a unit disk in the x-y plane to a square in 
plane with — 1 < [i < 1 and -1 < ( < 1. The simulation 
is then done in /z— £ plane with Cartesian coordinate. The 
resolution is mostly 256 x 256, with a few 512 x 512 runs 
to test convergence. By avoiding the polar coordinate 
and hence the singularity at the origin, we need a much 
smaller numerical viscosity term, iy\7 2 uj, to stabilize the 
simulation and hence obtain more reliable long time re- 
sults. It is also noted that, since it is impossible to use 
a true step vorticity distribution with finite grid points, 
the simulation results should not be compared exactly 
with the predictions based on Eq. (||). 

So for an initial ellipse with particular values of rg and 
e, we run simulations to long time and determine their fi- 
nal states. The results are plotted in Figure 1 as symbols, 
where squares indicate relaxation to elliptical states and 
solid circles to symmetrical states. Boundary region be- 
tween squares and circles represents the conditions where 
it is difficult to determine final states from simulations. 
In the figure we see the confirmation of stable elliptical 
vortices with large vortex sizes. All squares to the right of 
the solid line shows that simulations are consistent with 
the predictions. The simulations also show that ellipses 
are actually stable to a lower radius and the smallest 
stable size is decreasing with increasing e. Finally no 
relaxations to off-center vortices ever happen. 

In conclusion we have showed that, from vorticity mix- 
ing in time evolution and energy calculations, large ellip- 
tical vortices in a finite disk will remain stable. At the in- 
finitesimal ellipticity limit, this indicates stable I = 2 dio- 
cotron modes for large vortices. Numerical simulations 



not only confirmed these results, but also shows that el- 
liptical states are actually stable to a smaller size. The 
contradiction to current general idea of decaying I = 2 
modes is also indicated due to the incompleteness for con- 
sidering only monotonic decreasing vorticity by Briggs, 
Daugherty, and Levy. 

The author thanks Dr. C. Y. Lu for fruitful discussions 
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support of National Science Counsel, Taiwan, through 
the contract No. NSC 87-2112-M-008-034 is also ac- 
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